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\  Abstract 

The  Statistical  Package  for  the  Social  Sciences  (SPSS)  Discrimi¬ 
nant  Analysis  Routines  were  applied  to  speech  data  obtained  at  various 
gravity  levels.  The  data  were  created  on  the  Air  Force  Medical 
Research  Lab's  Centrifuge  at  Wright-Patterson  Air  Force  Base.  They 
were  then  digitally  sampled  and  Fourier  transformed  to  64  frequency 
bands  which  were  converted  to  two  separate  files;  one  the  16  frequencies 
used  by  Guyote  and  Sisson;  one  the  19  frequencies  used  by  Threshold 
Technology,  Inc.  These  files  were  then  further  processed:  summing 
the  energy  across  time  slices,  sutrming  across  frequency  bands,  and 
averaging  across  frequency  bands.  Statistical  results  showed  that 
G-stress  seemed  to  be  moving  energy  across  time  slices,  more  than 
shifting  across  frequencies. 
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APPLICATION  OF  DISCRIMINANT  ANALYSIS  TECHNIQUES 
TO  G-STRESS  ANALYSIS 
IN  COMPUTER  VOICE  DECODING 


I .  Introduction 

Background 

Despite  the  best  efforts  of  psychologists  and  aviators  to  lay  out 
control  panels  in  an  easy-to-read  display,  the  amount  of  visual  infor¬ 
mation  available  to  a  pilot  is  overwhelming,  particularly  in  a  one 
person  fighter  aircraft  where  the  pilot  may  need  to  do  more  than  keep 
the  plane  level.  In  addition,  the  number  of  buttons,  knobs,  levers 
and  switches  necessary  to  control  an  aircraft  keeps  the  plane  from 
reacting  at  any  speed  greater  than  the  manual  dexterity  of  the  pilot. 
With  the  introduction  of  a  new  channel  of  information  transfer,  i.e., 
speech,  the  pilot  could  both  process  more  information  and  control  more 
aircraft  functions  at  the  same  time,  thereby  increasing  the  plane's 
overall  reaction  capabilities. 

Although  the  field  of  speech  recognition  has  advanced  to  the 
point  where  isolated  words  can  be  recognized  with  a  99  percent  accuracy 
under  carefully  controlled  conditions  (Ref  2),  the  aircraft  environment 
has  many  speech  distortion  factors,  among  them:  G-forces,  oxygen  mask, 
noise,  vibration,  and  pilot  stress  and  fatigue. 

Threshold  Technology  has  a  voice  recognition  system  with  99.32 


percent  recognition  rate  for  normal  voiced  words  and  a  better  than  97 
percent  accuracy  rate  for  commands  spoken  through  a  standard  A.F. 
oxygen  mask  and  microphone  (Ref  10).  However,  there  are  no  systems  yet 


capable  of  recognizing  voiced  speech  distorted  by  high  G-forces  with 
any  reliable  degree  of  accuracy. 


Research  done  at  Rome  Air  Development  Center  from  5  January  1976 
to  4  January  1977  using  data  gathered  in  the  Brooks  AFB  centrifuge 
verified  the  existence  of  the  6-stress  problem,  but  failed  to  find 
any  pattern  in  the  distortion  (Ref  4). 

Objective 

The  objective  of  this  project  was  to  apply  the  Statistical  Pack¬ 
age  for  the  Social  Sciences  (SPSS)  Discriminant  Analysis  Techniques 
(Ref  6)  to  data  obtained  under  G-stress,  and,  using  the  results,  to 
form  a  pattern  of  the  G-stress  distortion. 

Basic  Approach 

To  use  the  SPSS  package  on  speech  data  an  immediate  problem  had 
to  be  surmounted.  Speech  data  are  often  processed  in  spectrographic 
form,  three  dimensional  arrays  which  can  be  pictured  as  2-D  graphs  of 
frequency  versus  time  with  intensity  as  the  third  dimension.  SPSS 
Discriminant  Analysis  is  designed  to  work  on  linear  vectors.  Cutting 
the  spectrographs  into  series  of  vectors  would  not  solve  the  problem 
because  the  number  of  variables  would  be  too  great  for  the  number  of 
samples  of  data.  SPSS  would  simply  truncate  the  data  arbitrarily. 

To  solve  this  problem,  the  spectrographic  data  were  "linearized" 
in  the  following  way  (see  Figure  1).  The  amplitudes  of  all  frequency 
components  within  a  time  slice  were  added  together  forming  a  vector  of 
Sum  of  Time  Slice  (SUMT).  Likewise  the  amplitude  of  all  the  time 
slices  for  a  single  frequency  were  summed  yielding  a  Sum  of  Frequency 


Vector  (FSUM).  Since  every  word  had  its  own  length,  the  FSUM  vectors 
were  partially  influenced  by  time,  and  so  a  third  vector  was  created, 
the  Average  of  Frequency  Vector  (FAVE),  by  dividing  each  FSUM  vector 
by  the  number  of  time  slices  over  which  it  was  summed.  Each  of  these 
vector  types  became  the  basic  record  for  a  file.  The  Discriminant 
Analysis  Tests  were  run  against  each  of  these  files  and  the  results 
compared  to  find  the  G-stress  effect. 


Figure  1.  Linearization  of  Spectrographic  Data 


The  discarding  of  information  in  each  of  these  vector  files 
causes  a  type  of  information  "hiding"  which  is  the  key  to  the  results 
of  this  thesis.  It  is  best  explained  by  imagining  that  all  of  the 
information  contained  in  the  speech  data  could  be  classified  into  say 
ten  distinct  components.  Now  suppose  ten  files  are  created,  each  one 
missing  one  of  those  components  of  information,  and  Discriminant 
Analysis  tests  are  run  for  each  of  the  files.  If  the  results  of  those 
tests  showed  that  one  file  had  no  distinctions  between  the  data  at 
different  gravity  levels  it  would  clearly  indicate  that  the  G-stress 
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was  affecting  only  that  one  component  of  information  that  was  missing. 
In  the  data  files  used  in  this  thesis,  the  FSUM  files  are  missing 
information  about  the  relative  timing  of  energy  amplitudes.  The  FAVE 
files  are  missing  this  and  the  total  energy/length  of  word  information 
The  SUMT  files  are  missing  information  on  the  relative  frequency 
amplitude  ratios. 

Scope 

The  scope  of  this  project  was  limited  to  seven  single  word 
comnands  spoken  by  two  volunteers  at  three  basic  gravity  levels.  The 
seven  words  used  in  this  project  were: 

WAITING,  BENCH,  FATHERS,  MY,  OUT,  AT,  WAS 
The  gravity  levels  were  IG,  1.5G  and  5G  for  the  first  speaker  and  IG, 
1.5G  and  4G  for  the  second  speaker. 

The  data  were  broken  into  two  data  sets.  The  first  included  only 
the  first  three  words  (WAITING,  BENCH,  FATHERS)  spoken  by  the  first 
speaker  (IG,  1.5G,  5G).  The  Discriminant  Analysis  results  obtained 
from  this  limited  data  suggested  the  pattern  of  energy  change  over 
time.  The  second  data  set  was  used  to  further  test  this  theory.  It 
included  the  first  three  words  spoken  by  the  second  speaker,  the  last 
four  words  spoken  by  the  first  speaker,  and  all  seven  words  spoken 
by  a  third  speaker  not  wearing  an  oxygen  mask.  The  data  sets  seem 
small,  however,  each  one  completely  filled  a  7-track,  2400  ft,  800 
bytes  per  inch  tape  when  digitized. 

All  data  input  from  the  first  two  speakers  were  obtained  using 


the  standard  A.F.  oxygen  mask  and  microphone.  Although  the  mask 
itself  creates  a  speech  distortion,  the  presence  of  the  mask  at  all 


II.  Data  Acauisition 


Data  acquisition  was  performed  in  the  same  basic  manner  as 
described  by  Guyote  and  Sisson  (Ref  1). 

Raw  Data  Acquisition 

The  original  data  tape  used  in  this  research  was  supplied  by 
Mr.  Don  McKecknie  of  the  Air  Force  Medical  Research  Lab  of  Wright- 
Patterson  Air  Force  Base.  Two  volunteers  were  recorded  in  the  centri¬ 
fuge  on  two  different  channels  of  the  same  tape.  The  quality  of  this 
tape  was  somewhat  less  than  ideal  because  of  a  high  degree  of  breath 
noise  in  the  oxygen  mask  microphone. 

Each  of  the  speakers  recorded  the  same  set  of  twelve  words: 

WAITING,  BENCH,  FATHERS,  LAWN,  SHE, 

TOOK,  SHOE,  JOE,  MY,  OUT,  AT,  WAS 

They  were  recorded  first  at  IG,  with  the  centrifuge  standing  still, 
then  at  1.5G,  considered  a  baseline  level  for  the  centrifuge.  Each 
volunteer  was  then  given  control  for  the  high  gravity  run.  The  first 
speaker  took  the  centrifuge  to  5G  and  was  able  to  repeat  two  samples 
of  each  word  duri*"  his  30  second  run.  The  second  speaker  recorded 
three  samples  each  at  4G. 

This  original  tape  was  then  edited  with  the  use  of  an  entertain¬ 
ment  grade  tape  cassette  player  and  re-recorded  on  an  unused  portion 
of  the  original  1/4  inch  tape.  The  use  of  this  equipment  was 
justified  since  the  end  result  sounded  the  same  to  the  human  ear, 
indicating  that  all  the  information  necessary  for  word  recognition 
remained.  In  fact,  the  words  used  in  the  limited  data  space  available 
were  determined  by  their  clarity  to  the  human  ear.  For  ease  of 
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recognition  of  critical  points  in  the  tape,  a  short  harmonica  tone  was 
recorded  between  speakers  and  gravity  levels.  These  notes  showed  up 
very  clearly  on  the  spectrogram  isolating  data  words  from  the  surround¬ 
ing  noise. 

Analoq-to-Diqital  Conversion  and  FFT 

The  edited  data  were  digitized  by  Mr.  Jack  Capehart  of  ASD/AD 
following  the  same  procedures  originated  by  Maj  Neyman  (Ref  5),  These 
orocedures  resulted  in  an  effective  sampling  rate  of  10  kHz  with  the 
speech  signal  low  pass  filtered  to  5  kHz. 

The  digitized  speech  was  then  Fast  Fourier  Transformed  (FFT) 
using  a  window  of  12.8  ms,  or  128  samples  per  FFT,  This  resulted  in 
a  set  of  64  discrete  amplitude  values,  representing  64  frequency 
multiples  of  78.125  Hz  for  each  12.8  ms  time  slice.  Three  minutes 
worth  of  this  data  could  be  stored  on  one  tape. 

Spectrogram  and  Data  Selection 

Spectrographic  printout  of  the  entire  data  tape  was  accomplished 
using  a  Fortran  program  provided  by  Capt  William  Nelson,  AFIT/EN 
(see  Appendix).  This  program  was  an  improvement  over  the  one  used  by 
Maj  Neyman  (Ref  5)  in  that  it  used  only  a  single  overprint  to  create 
the  spectrographic  light/dark  effect.  This  program  was  also  flexible 
enough  to  be  used,  with  only  minor  variations,  on  the  various  data 
files  created  from  the  original  tape. 

The  data  representing  the  spoken  words  were  then  isolated  from 


the  surrounding  noise.  Beginning  and  end  points  of  words  were  deter¬ 
mined  by  visual  study  of  the  data  tape  spectrogram  (see  Figure  2) 
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Figure  2.  Sample  Spectrograph 
(Word  is  'MY'  at  IG  with  no  mask.  Data  Set  2) 
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which  included  tape  record  numbers  and  sum  of  squares  of  the  frequency 
amplitudes  for  each  record.  Precedence  for  this  procedure  can  be  found 
in  Reference  7. 

Program  GETDATA  (see  Appendix)  then  created  the  first  main  data 
file  (named  TALKDATA),  stored  on  disk  for  fast  access.  It  did  this 
by  alternately  bypassing  groups  of  records  on  the  tape  representing 
noise,  and  copying  groups  of  records  from  the  tape  which  contained 
the  spoken  word  data.  (Each  record  on  the  tape  contained  64  Fourier 
Transform  values  representing  one  12.8  ms  time  slice.)  Within  the 
copying  routine,  a  call  to  subroutine  MAKEPF  "standardized"  the  size 
of  each  word.  For  example,  suppose  one  utterance  of  a  word  at  a 
given  gravity  level  lasted  28  time  slices  and  another  utterance 
lasted  33  time  slices.  The  first  word  would  have  two  time  slices 
created  by  averaging  time  slices  10  and  11  and  20  and  21.  The  second 
would  have  time  slices  9,  17  and  25  removed.  Both  words  would  then 
be  exactly  30  time  slices  long,  the  template  size  representing  that 
word  and  gravity. 
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III.  Data  Base  Creation 


Data  Files  and  Structure 

From  the  64  channel  file  TALKDTA,  two  smaller  files  were  created: 
TALK16  with  16  frequency  channels,  and  TALK19  with  19  frequency 
channels  (see  Table  I).  These  files  were  created  both  because  64 
channels  contained  too  many  variables  for  the  SPSS  package  to  handle 
properly  with  the  available  number  of  data  samples,  and  because  the 
amplitudes  of  the  higher  frequencies  were  so  small  compared  to  the 
amplitudes  of  the  lower  frequencies,  that  they  could  not  be  effectively 
used  as  discriminating  variables. 

TALK16  was  created  by  program  CUTDOWN  (see  Appendix)  using  the 
same  methods  as  Guyote  and  Sisson  (Ref  1).  It  logarithmically 
increases  the  values  of  the  higher  frequency  amplitudes  while  cutting 
down  the  number  of  frequency  channels.  This  set  of  16  frequency 
channels  was  selected  for  study  in  this  project  because  of  its  common 
use  in  other  speech  recognition  schemes. 

TALK19  was  created  by  program  CONVERT  (see  Appendix).  This  set 
of  19  frequency  channels  was  selected  for  study  because  Threshold 
Technology  declared  it  to  be  the  best  set  of  input  frequencies  for 
its  speech  recognizer,  which  has  been  used  successfully  in  speech 
recognition  tests  where  the  speakers  were  wearing  oxygen  masks  (Ref  8). 

The  following  procedures  were  used  to  create  file  TALK19: 

1.  The  Threshold  Technology  Voice  Control  Demonstration 
System  Operating  and  Maintenance  Manual  (Ref  9)  displayed  energy  out¬ 
puts  for  their  filters  at  a  given  frequency.  From  this  information,  a 
curve  representing  the  frequency  response  of  each  frequency  could  be 
plotted. 


Table  I 


19  Channels 

16  Channels  (used  by  Threshold  Technology  Voice 

(as  per  1977  Thesis  by  Control  Demonstration  System  Equipments 

M.  F.  Guyote  and  P.  L.  Sisson)  Manual  Schematic  &  Parts  List 


1. 

78.125 

Hertz 

1. 

260 

Hertz 

2. 

156.250 

Hertz 

2. 

395 

Hertz 

3. 

234.375 

Hertz 

3. 

535 

Hertz 

4. 

312.500 

Hertz 

4. 

683 

Hertz 

5. 

390.625 

Hertz 

5. 

841 

Hertz 

6. 

468.750 

Hertz 

6. 

1011 

Hertz 

7. 

585.940 

Hertz 

7. 

1198 

Hertz 

8. 

742.188 

Hertz 

8. 

1405 

Hertz 

9. 

898.440 

Hertz 

9. 

1635 

Hertz 

10. 

1132.810 

Hertz 

10. 

1892 

Hertz 

11. 

1445.310 

Hertz 

11. 

2179 

Hertz 

12. 

1793.380 

Hertz 

12. 

2505 

Hertz 

13. 

2226.560 

Hertz 

13. 

2885 

Hertz 

14. 

2812.500 

Hertz 

14. 

3326 

Hertz 

15. 

3554.690 

Hertz 

15. 

3855 

Hertz 

16. 

4453.125 

Hertz 

16. 

4484 

Hertz 

17. 

5263 

Hertz 

18. 

6277 

Hertz 

19. 

7626 

Hertz 

2.  Using  standard  numeric  analysis  curve  fitting  techniques, 
the  relative  amplitudes  (from  0  to  1)  of  the  basic  64  frequency 
channels  could  be  calculated  and  summed. 

3.  The  resulting  19  amplitudes  represent  what  the  Threshold 
Technology  filters  would  have  recorded  given  the  same  speech  input 
represented  by  the  64  channel  file. 

Program  SUMFREQ  (see  Appendix)  was  then  used,  once  for  each  file 
TALK16  and  TALK19,  to  create  the  final  Data  Base  Sets.  Word  parameters 
were  read  in  first  (WORD,  GRAVITY,  SPEAKER,  and  WORD  LENGTH)  followed 
by  the  time/frequency  amplitudes  in  their  Number  of  Frequency  x  Number 
of  Time  Slices  data  arrays.  These  arrays  were  then  summed  and 
averaged  as  explained  in  the  introduction  (see  Figure  1  again).  The 
final  resulting  files  were  labeled  FSUM16,  FAVE16,  SUMT16,  FSUM19, 
FAVE19,  SUMT19  for  the  first  data  set.  The  second  data  set  consisted 
of  FSUM216,  FAVE216,  SUMT216,  FSUM219,  FAVE219,  SUMT219  (the  "2" 
signifying  second  data  set). 

In  the  first  data  set,  the  first  and  third  words  were  30  time 
slices  long.  The  second  word  was  only  18  time  slices  long,  but  was 
zero  filled  to  create  30  numbers  for  the  SUMT  vectors.  In  the  second 
data  set,  the  longest  word  was  60  time  slices  and  so  60  became  the 
template  size.  All  other  SUMT  data  words,  which  had  shorter 
"standard"  time  slice  duration,  were  fitted  into  the  beginning  of 
the  60  time  slice  template,  with  zeros  used  for  fillers. 


IV.  Discrimi nant  Analysis 

Discriminant  Analysis  Capabilities 

Discriminant  analysis  is  a  statistical  method  of  determining 
distinctive  group  characteristics  of  known  grouped  data,  and  deter¬ 
mining  group  membership  of  ungrouped  data  according  to  those  character¬ 
istics.  The  way  in  which  this  is  done  is  analogous  to  finding  a  view 
of  all  of  the  data  which  best  separates  the  groups.  Figure  3a 
represents  one  dimensional  data  such  as  age  of  employees,  or  salary 
levels.  ...  In  the  case  pictured,  the  dimension  can  easily 
discriminate  between  the  two  groups  of  data,  those  above  zero  and 
those  below  zero.  Any  new  data  can  then  be  grouped  as  positive  or 
negative. 

Figure  3b  represents  data  consisting  of  two  variables,  conven¬ 
iently  labeled  x  and  y.  Notice  that  if  only  x  or  y  were  used  to 
discriminate  among  the  groups  there  would  be  much  group  overlapping 
and  group  membership  classification  would  be  difficult.  However,  if 
X  and  y  are  used,  the  group  memberships  become  clearly  distinct. 

Figure  3c  represents  data  sets  of  three  dimensions.  In  the  case 
pictured,  neither  x  nor  y  can  help  discriminate  between  the  groups, 
but  the  z  data  variable  alone  can  clearly  define  group  membership. 

Clearly,  data  sets  can  have  any  "n"  number  of  variables,  and  may 
be  split  into  any  "g"  number  of  groups.  The  Discriminant  Analysis 
Program  can  determine  which  variable  best  distinguishes  among  the 
groups  and  how  much  distinction  it  makes.  It  can  then  bring  in  the 
next  most  useful  variable  for  group  classification,  and  so  on  until 
either  all  the  variables  are  used,  the  number  of  variables  remaining 
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add  nothing  to  the  discrimination,  or  the  number  of  variables  allowed 
has  been  reached. 

There  are  many  different  ways  in  which  the  distinguishing  vari¬ 
ables  can  be  chosen.  In  this  thesis  each  test  was  run  using  two 
different  methods:  MAHAL  and  MAXMINF.  This  helped  to  determine  if 
one  method  of  selection  was  better  than  the  other  for  any  purpose, 
or  whether  the  same  variables  would  be  consistently  chosen  as  most 
discriminating.  Method  =  MAHAL  chooses  variables  that  maximize  the 
mahalanobis  distance  (similar  to  a  Euclidian  sum  of  squares)  between 
group  centroids.  Method  =  MAXMINF  chooses  variables  that  maximize 
MINF,  where  F  is  a  standard  statistical  two  group  F  function  and 
MINF  is  the  minimum  F  value  among  all  the  inter-group  F  values.  (A 
standard  two  group  F  function  measures  the  ratio  of  the  variance  of 
two  groups.  Variance  is  a  measure  of  the  spreading  or  grouping  of 
data  on  a  graph  such  as  in  Figures  3a,  3b  and  3c  which  show  tightly 
grouped  data.) 

Tests  Selected 

The  significance  of  tests  was  originally  determined  through 
trial  and  error.  Initial  efforts  were  designed  to  compare  each  word 
to  itself  at  differing  gravity  levels,  and  then  compare  the  resultant 
discriminant  and  classification  functions  to  each  other.  However,  no 
pattern  was  discerned  from  this  method. 

Next,  the  different  words  were  compared  against  each  other  at 
each  of  the  same  gravity  levels.  The  words  at  different  gravity 
levels  were  then  grouped  according  to  the  results  of  the  classification 
functions  derived  from  the  determining  set.  From  these  tests  it  was 
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noted  that  the  results  for  file  FSLIM16  were  particularly  good;  good 
in  that  the  classification  functions  derived  from  any  one  gravity 
level  set  of  data  were  able  to  correctly  identify  the  words  of  the 
other  gravity  levels. 

The  tests  were  then  refined  to  limit  the  number  of  variables 
allowed  for  use  in  the  classification  functions.  This  was  done  both 
to  increase  the  validity  of  the  discrimination,  and  to  determine 
which  variables,  if  any,  would  be  chosen  consistently  as  discriminating 
characteristics. 

Tests  were  also  conducted  to  measure  the  distinguishability  of 
each  word  if  all  gravity  levels  were  used  as  a  basis;  and  to  separate 
each  individual  word-gravity  combination  as  much  as  possible. 

Test  Specifications  (Data  Base  1 ) 

Before  the  Group  Discrimination  Tests  could  be  run,  the  data 
first  had  to  be  separated  into  groups.  One  obvious  grouping  simply 
used  the  variable  IWORD  to  determine  if  a  data  record  represented 
Word  1,  2  or  3.  Another  possible  grouping  would  have  been  to  divide 
the  records  (using  IGRAV)  into  three  groups  based  on  gravity  level 
alone. 

Since  it  was  possible  to  compute  new  functions  for  the  purpose 
of  group  identification,  the  data  records  were  also  separated  into 
nine  groups  representing  each  unique  word/gravity  pair  (3  words  x 
3  gravity  levels  =  9  unique  possible  groups  for  Data  Set  1).  Two 
different  naming  functions  were  used  to  make  reading  of  the  SPSS 
results  easier. 


The  function  IGRAVWORD  =  (IGRAV  -  1)  x  3  +  IWORD  sets  up  nine 
groups  which  are  semi -grouped  by  gravity  level.  They  are: 

Groups  1,  2,  3  =  Words  1,  2  &  3  at  IG 

Groups  4,  5,  6  =  Words  1,  2  &  3  at  1.5G 

Groups  7,  8,  9  =  Words  1 ,  2  &  3  at  5G 

The  function  IWORDGRAV  =  (IWORD  -  1)  x  3  +  IGRAV  sets  up  groups 
which  are  semi-grouped  by  word  so  that: 

Groups  1,  2,  3  =  Word  1  at  1,  1.5  and  5G 

Groups  4,  5,  6  =  Word  2  at  1,  1.5  and  5G 

Groups  7,  8,  9  =  Word  3  at  1 ,  1.5  and  5G 

The  test  IWORDGRAV  (1,9)  set  the  Discriminant  Analysis  Program 
to  separate  out  each  of  the  nine  words/gravities  as  distinctly  as 
possible.  If  the  resultant  groupings  left  the  same  word  at  different 
gravity  levels  close  together,  it  would  tend  to  indicate  that  the 
G-stress  had  little  effect  on  that  data  set.  Such  a  result  should 
be  accompanied  by  a  low  error  rate  on  the  IWORD  (1,3)  test,  which 
classifies  the  3  word  groups  using  only  word  classification  for 
discrimination.  Close  grouping  of  the  words  at  different  gravity 
levels  would  lead  to  a  low  error  rate,  and  would  be  due  to  a  low 
G-stress  effect. 

The  test  IGRAVVIORD  (1,3)  compares  IGRAVWORD  groups  1,  2  and  3 
to  each  other  and  produces  a  set  of  classification  functions  which 
can  discriminate  among  the  three  words  at  gravity  level  1.  SPSS  then 
classifies  all  the  other  word  data  sets  according  to  these  same 
classification  functions.  These  are  called  "ungrouped"  data  sets 
because  they  do  not  belong  to  groups  1,  2  or  3.  Tests  IGRAVWORD  (4,6) 


and  IGRAVWORD  (7,9)  perform  the  same  function  as  IGRAVWORD  (1,3) 
using  the  gravity  levels  1.5G  and  5G  respectively  as  the  basis  for 
classification. 
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V.  Discriminant  Analysis  Initial  Results  (Data  Base  1) 

The  results  of  the  Discriminant  Analysis  Tests  on  Data  Base  1  are 
summarized  in  Table  II.  The  term  VAR#####  is  a  list  of  the  variables 
used  in  the  classification  function.  The  term  ERR  refers  to  a  word 
being  incorrectly  classified  in  any  way.  The  term  WW  stands  for  Wrong 
Word  and  indicates  that  some  word  was  misclassified  as  another.  It  is 
listed  as  a  fraction  in  Part  II  of  Table  II  indicating  the  number 
of  words  classified  incorrectly/total  number  of  ungrouped  words 
classified. 

Note  that  an  error  can  occur  that  is  not  a  wrong  word  if,  for 
example,  a  word  at  IG  is  classified  as  being  the  same  word  at  5G.  The 
implication  of  this  happening  is  that  the  discrimination  program  could 
not  separate  a  word  from  itself  at  a  different  gravity  level  although 
it  could  easily  tell  which  word  it  was.  Alternately,  a  wrong  word 
indicates  that  the  gravity  so  changed  a  word  that  it  was  unrecognizable 
as  itself.  The  effect  of  the  G-stress  on  each  of  the  Data  Base  1 
files  could  thus  be  measured  as  a  function  of  the  number  of  errors 
and  the  number  of  wrong  words. 

As  can  be  seen,  the  results  of  test  IWORDGRAV  (1,9)  indicate  a 
distinctly  higher  error  rate  for  the  files  SUMT16  and  SUMT19  than  for 
the  other  four  files.  The  results  of  the  Part  II  tests  are  not  as 
marked.  The  difference  in  error  rates  between  FAVE16,  FSUM16,  FSUM19 
and  SUMT16,  SUMT19  follows  the  same  pattern  as  the  Part  I  tests, 
however,  the  error  rate  for  FAVE19  is  too  high  to  fit  the  pattern. 
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table  U 

Results  (continued) 


VI.  Data  Base  1  Conclusions 

The  high  number  of  wrong  words  selected  for  files  SUMT16  and  SUMT19 
in  Part  I  indicates  that  the  gravity  word  distinction  affect  is  stronger 
in  these  two  files  than  in  the  other  four  files.  This  is  exactly  the 
opposite  of  what  would  happen  if  gravity  was  only  shifting  energy  from 
one  frequency  to  another.  Therefore,  although  the  frequency  shift  is 
not  ruled  out,  it  seems  that  energy  is  shifting  across  time  slices  to 
a  greater  extent  than  any  frequency  shifting  that  might  be  present. 

And  even  though  the  number  of  errors  which  are  not  wrong  word  choices 
is  also  high  (which  would  indicate  a  similarity  across  the  gravity 
levels),  the  Wrong  Word  rate  for  SUMT16  and  SUMT19  is  50  percent  of 
all  errors  whereas  for  the  three  files  FAVE16,  FSUM16  and  FSUM19  the 
Wrong  Word  rate  is  less  than  31  percent. 

The  same  conclusions  can  be  reached  from  the  error  rate  results  of 
Part  II.  The  low  error  rates  for  files  FAVE16,  FSUM16  and  FSL)M19 
indicate  that  classifications  of  words  grouped  at  any  one  gravity  will 
hold  for  those  words  at  different  gravities,  indicating  small  gravity 
distortion  across  the  frequency  bands.  The  high  error  rate  of  SUMT16 
and  SUMT19  indicate  gravity  distortion  across  time. 

A  possible  explanation  for  this  effect  might  be  that  at  higher 
gravity  levels,  the  pilot  must  force  speech  out  in  a  sharp  exhalation 
of  breath.  This  might  cause  greater  energy  emphasis  at  the  beginning 
of  an  utterance,  or  otherwise  alter  the  timing  of  the  components  of 
the  utterance. 

File  FAVE19  is  an  enigma.  In  Part  I  it  has  a  wrong  Word  ratio  of 
100  percent,  and  in  Part  II  it  has  a  high  error  rate,  behaving  just 
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like  the  files  SUMT16  and  SUMT19.  But  since  its  data  is  completely 
time  independent  (representing  the  average  energy  per  time  slice  of 
each  frequency)  this  behavior  does  not  fit  the  pattern  of  the  other 
five  files,  and  challenges  the  conclusions  based  on  that  pattern. 

It  was  hoped  that  the  results  of  the  second  data  set  would  be 
according  to  the  first  pattern.  But  if  the  FAVE219  data  behaved  as 
file  FAVE19  did,  then  the  theory  of  G-stress  distortion  would  have 
to  take  this  consistency  into  account  and  determine  why  this  file 
behaved  thus. 


23 


VII .  Data  Base  2 


Data  Structure 

Data  files  FAVE216,  FSUM216,  SUMT216,  FAVE219,  FSUM219  and  SUMT219 
were  created  in  a  similar  manner  to  Data  Base  1.  The  main  difference, 
physically,  is  that  some  of  the  word  data  in  this  data  base  extended 
for  60  time  slices,  which  when  summed  within  time  slices  led  to  a 
SUMT216  and  SUMT219  record  greater  than  150  characters.  Since  SPSS 
cannot  handle  such  large  records,  these  two  files  had  to  be  rewritten 
as  two  record  data  sets.  It  was  also  necessary  to  sort  the  second 
data  set  in  order  to  more  easily  interpret  the  SPSS  results.  Neither 
of  these  file  differences  had  any  effect  on  the  data. 

The  main  difference  between  the  two  data  sets,  logically,  is  that 
Data  Set  II  contains  seven  words,  three  speakers  and  five  gravity 
levels  (gravity  level  5  is  actually  IG  with  the  speaker  wearing  no 
mask). 

Following  is  a  list  of  Data  Base  2: 

Word  1  =  WAITING  Speaker  2  at  G-levels  1,  2,  4 
and  Speaker  3  at  G-level  5 

Word  2  =  BENCH  Speaker  2  at  G-levels  1,  2,  4 
and  Speaker  3  at  G-level  5 

Word  3  =  FATHERS  Speaker  2  at  G-levels  1,  2,  4 
and  Speaker  3  at  G-level  5 

Word  4  =  MY  Speaker  1  at  G-levels  1,  2 

and  Speaker  3  at  G-level  5 

Word  5  =  OUT  Speaker  1  at  G-levels  1,  2,  3 

and  Speaker  3  at  G-level  5 

Word  6  =  AT  Speaker  1  at  G-levels  1,  2,  3 

and  Speaker  3  at  G-level  5 

Word  7  =  WAS  Speaker  1  at  G-levels  1,  2,  3 

and  Speaker  3  at  G-level  5 
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Tests  Selected  and  Specifications 


Function  IGRAVWORD  =  [(IGRAV  -  1)  *  7  +  IWORD]  sets  up  the  35 


groups ; 

Groups  1-7  =  all  seven  words  at  Gravity  Level  1  (IG) 

Groups  8-14  =  all  seven  words  at  Gravity  Level  2  (1.5G) 

Groups  15-21  =  all  seven  words  at  Gravity  Level  3  (5G) 

Groups  22-28  =  all  seven  words  at  Gravity  Level  4  (4G) 

Groups  29-35  =  all  seven  words  at  Gravity  Level  5  (IG  -  no 

mask) 

Function  IWORDGRAV  =  (IWORD  -  1)  *  5  +  IGRAV  sets  up  the  35 


groups: 

Groups  1-5  =  Word  1  at  all  Gravity  Levels 

Groups  6-10  =  Word  2  at  all  Gravity  Levels 

Groups  11-15  =  Word  3  at  all  Gravity  Levels 

Groups  16-20  =  Word  4  at  all  Gravity  Levels 

Groups  21-25  =  Word  5  at  all  Gravity  Levels 

Groups  26-30  =  Word  6  at  all  Gravity  Levels 

Groups  31-35  =  Word  7  at  all  Gravity  Levels 

Note  that  not  every  group  actually  has  data  in  it.  For  instance, 
there  is  no  data  for  IGRAVWORD  groups  25-28  because  speaker  2  only 
spoke  Words  1,  2  and  3  and  no  other  speaker  was  recorded  at  4G's. 
This  does  not  invalidate  the  test  results.  Whenever  a  group  has  no 
data,  SPSS  simply  leaves  it  out  of  the  classification  function. 

The  tests  run  on  Data  Base  2  were  selected  to  parallel  the  Data 
Base  1  tests.  Test  IWORDGRAV  (1,35)  is  the  same  as  Data  Set  I's 
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IW0RD6RAV  (1,9)  in  that  it  discrininates  each  word/gravity  group  from 
the  other.  The  much  higher  error  rate  of  this  test  in  the  second  data 
base  is  not  because  of  the  greater  number  of  words,  but  rather  the 
fact  that  each  word  was  spoken  by  two  speakers;  one  with  no  face  mask 
at  IG,  and  the  other  with  a  face  mask  at  various  G-levels.  The  high 
error  rate  indicates  that  these  words  together  do  not  form  a  coherent 
group  distinctly  different  from  the  other  word  groups. 

In  contrast,  the  tests  IGRAVWORD  (x,y)  all  compare  all  seven 
words  spoken  at  the  same  gravity  level  (and  the  same  speaker).  This 
compares  exactly  to  the  Data  Base  1  tests. 
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IX.  Discriminant  Analysis  Results  (Data  Base  2) 

Table  III  shows  the  results  of  SPSS  testing  on  Data  Base  2.  As 
before,  WW  stands  for  wrong  word,  ERR  for  errors,  and  VAR  for  variables 
chosen,  and  the  results  are  in  fractions,  number  wrong/number  available. 

An  extra  level  of  complexity  has  been  added  to  the  Data  Base  2 
results  table  due  to  the  number  of  speakers.  That  is  why  there  are  two 
error  fractions  for  each  test.  Each  is  labeled,  either  same  speaker 
error  rate  or  other  speaker  error  rate  (at  the  same  gravity  level  or 
different  gravity  level).  The  totals  have  also  been  separated  out  in 
this  manner,  and  the  number  of  errors  in  the  grouped  data  have  also 
been  accounted  for. 

The  results  of  test  IWORDGRAV  (1,35)  (Part  I)  follow  the  same 
pattern  as  for  Data  Base  1.  Both  SUMT216  and  SUMT219  have  a  WW/ERR 
ratio  of  74  percent,  while  FSUM216  and  FSUM219  have  a  WW/ERR  ratio 
of  46  percent.  FAVE216  and  FAVE219  are  66  percent  and  63  percent 
respectively. 

Likewise,  the  results  of  Part  II  totals  for  the  same  speaker 
also  follow  the  same  basic  pattern  of  Data  Set  I,  with  SUMT216  and 
SUMT219  having  distinctly  higher  error  rates  than  FSUM216  and  FSUM219. 
However,  the  error  rates  for  FAVE216  and  FAVE219  are  only  marginally 
lower  than  the  error  rates  for  SUMT216  and  SUMT219.  It  can  also  be 
seen  that  the  error  rate  for  a  different  speaker  (with/without  fave 
mask)  is  approximately  twice  the  gravity  difference  error  rate  for  the 
same  speaker.  And  a  change  of  gravity  level  at  that  point  does  not 
raise  the  error  rates. 
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Table  lit  (continued) 


Data  Base  2  Conclusions 


The  results  of  Part  I  test  I  (discriminating  each  word/gravity 
pair  from  all  the  others)  and  the  results  of  Part  II  totaled  for  the 
same  speaker  tend  to  support  the  theory  that  gravity  is  affecting 
energy  across  time  slices  more  than  across  frequency  bands. 

However,  the  results  are  far  from  conclusive.  If  only  the  files 
FSUM216  and  FSUM219  vs  SUMT216  and  SUMT219  were  used  in  the  analysis 
the  results  would  be  clear.  However,  the  results  of  files  FAVE216 
and  FAVE219  do  not  behave  well  according  to  the  theory  proposed.  The 
possibility  exists  that  the  error  rates  for  FSUM216  and  FSUM219  are 
lower  than  the  other  files  because  they  retain  more  of  the  information 
available  in  the  original  word  spectrograms.  Certainly  files  FAVE216 
and  FAVE219  have  lost  the  information  of  the  word  length.  However, 
it  is  unclear  what  information  if  any  was  lost  to  files  SUMT216  and 
SUMT219  except  the  ability  to  discriminate  among  the  frequencies. 

Since  FSUM216  and  FSUM219  cannot  discriminate  across  time  slices,  it 
would  seem  that  their  relative  error  rates  would  indicate  a  valid 
pattern. 
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XI.  Recommendations 


Obviously,  more  data  are  needed  before  concrete  conclusions  can 
be  drawn.  The  quality  and  content  of  that  data  should  also  be  better 
controlled.  For  instance,  for  every  speaker  used  in  the  experiment, 
the  data  base  should  contain  speech  at  IG  -  no  mask,  speech  at  IG  - 
with  mask,  centrifuge  baseline  with  mask,  and  at  least  two  other 
gravity  level  rides  with  mask.  Also,  these  high  gravity  level  rides 
should  all  be  taken  at  the  same  gravity  level  for  all  riders. 

A  second  recommendation  is  that  more  research  be  done  on  other 
ways  to  set  up  the  data  files.  Specifically,  what  other  components 
of  the  speech  data  can  be  isolated  and  extracted  to  create  new  files 
with  new  components  of  information  hidden. 

A  third  recotmiendation  concerns  the  data  results  in  Tables  I,  II 
and  III.  The  variables  used  for  each  set  of  group  classifications  is 
listed  for  both  Data  Sets.  There  is  some  apparent  heavy  use  of 
frequency  5  for  FSUM216,  FAVE216,  FSUM16  and  FAVE16  and  a  corresponding 
(though  not  equal)  heavy  use  of  frequency  2  in  fi'--;  F5UM219,  FAVE219, 
FSUM19  and  FAVE19.  Channel  5  in  the  16  frequency  band  =  390.623  Hz. 
Channel  2  in  the  19  frequency  band  =  395  Hz.  This  shows  that  this 
frequency  is  a  main  discriminant  between  the  various  words.  Further 
study  might  reveal  other  key  frequencies.  These  frequencies  might 
then  form  a  basis  for  the  information  hiding  studies  just  discussed. 

Finally,  some  recommendations  for  work  that  could  be  performed 
on  the  data  bases  already  available  which  was  not  performed  simply 
because  of  a  lack  of  time.  Data  Sets  1  and  2  could  be  merged  into  a 
third  Data  Set  containing  the  first  3  words  spoken  at  all  available 
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gravity  levels  by  all  three  speakers.  This  would  make  it  possible  to 
do  preliminary  studies  on  the  effects  of  different  speakers  at  the 
same  gravity  level  with  both  wearing  oxygen  masks  as  compared  to  a 
different  speaker  without  a  face  mask.  (More  data  would  be  needed  for 
more  comprehensive  studies  of  the  separation  of  face  mask,  G-stress 
and  speaker  identity  effects.)  Also,  all  the  available  data  words 
could  be  stretched  or  cut  to  a  single  template  size  (no  filler)  using 
the  routines  in  program  GETDATA  (see  Appendix).  Not  only  could  furthe 
tests  be  conducted  on  this  data,  but  each  word  could  be  further 
normalized  to  have  the  same  overall  amount  of  energy  by  simply 
dividing  each  amplitude  in  each  template  by  the  energy  sum  total. 

This  might  eliminate  so  much  information  that  one  word  could  not  be 
distinguished  from  another,  or  it  might  lead  to  very  clearcut  word 
independent  G-stress  effects.  The  only  way  to  find  out  is  to  try  it. 
[Some  more  elegant,  but  very  time  consuming  computer  routines  have 
been  employed  by  (at  least)  NEC  and  Lear-Siegler  Corporation  to  handl 
this  problem  by  a  number  of  dynamic  programming  methods.  This  thesis 
would  use  a  simpler  approach  since  G-stress  is  the  prime  concern,  not 
word  recognition.] 
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PROGRAM  SUMFREQ  (INPUT, OUTPUT, TAPE1,TAPE2,TAPE3,TAPE4) 


CCC 

CCC  THIS  PROGRAM  WILL  CONVERT  FILE  TALK19  TO  FSUM19 
CCC  WHICH  WILL  TF  THE  SUM  OF  EACH  FREQUENCY  ACROSS 
CCC  TIME.  OUTPUT  FORMAT  WILL  BE:  IWORD,IGRAV, 

CCC  ISPKER,  IWRDSZ,  19  FREQ  SUMS. 

CCC  THIS  WILL  BE  SUITABLE  FOR  INPUT  TO  SPSS. 

CCC  ALSO  CREATED  WILL  BE  SUMT19,  WHICH  IS  THE 

CCC  SUM  OF  ALL  FREQUENCIES  WITHIN  A 

CCC  TIME  SLICE.  FORMAT  IS:  IWORD,IGRAV, ISPKER, 

CCC  IWRDSZ,  60  SUMT'S. 

CCC 

DIMENSION  A(60,19),  B(19),  C(60),  IC(60) 

CCC 

7.WORD  =  1 

ICRAV  =  1 

NBRTIM  =60 

DO  190  ISPEC  =1,90 

CCC 

CALL  READA  (A, NBRTIM, IWORD,IGRAV, ISPKER, IWRDSZ) 

CALL  SUM?  ( A, B, NBRTIM) 

CALL  SUMl  U.IC, NBRTIM) 

WRITE  (2,10)  IWORD,IGRAV, ISPKER, IWRDSZ, (B(J),  J  =  1,19) 
WRITE  (3,20)  IWORD,IGRAV, ISPKER, Il^RDSZ,(IC(I), 1=1, 60) 

WRITE  (4,10)  IWORD,IGRAV, ISPKER, IWRDSZ, ((B(J)/IWRDS2),J=1, 19) 
10  FORMAT  (4I1,19F6.1) 

20  FORMAT  (411,6014) 

CCC 

CCC  IF  (ISPEC.EQ.4)  IGRAV  =  2 
CCC  IF  (ISPEC.EQ.8)  IGRAV  =  3 
CCC  IF  ( ISPEC. EQ. 10)  IGRAV  =  1 

CCC  IF  ( ISPEC. EQ. 10)  IWORD  =  2 

CCC  IF  (ISPEC. EQ. 10)  NBRTIM  =  18 
CCC  IF  (ISPEC. EQ. 14)  IGRAV  =  2 

CCC  IF  ( ISPEC. EQ. 18)  IGRAV  =  3 

CCC  IF  ( ISPEC. EQ. 20)  IGRAV  =  1 

CCC  IF  ( ISPEC. EQ. 20)  IWORD  =  3 

CCC  IF  ( ISPEC. EQ. 20)  NBRTIM  =  30 
CCC  IF  ( ISPEC. EQ. 24)  IGRAV  =  2 

CCC  IF  ( ISPEC. EQ. 28)  IGRAV  =  3 

190  CONTINUE 
STOP 
END 

CCC 

CCC 
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SUBROUTINE  READA  (A, NBRTIM, IWORD , IGRAV, ISPKER, IWRDSZ) 
DIMENSION  A(60,19) 

READ  (1,*)  IWORD, IGRAV, ISPKER, IWRDSZ 
DO  100  I  =  1, NBRTIM 
DO  100  J  =  1,19 
READ  (1,*)  A(I,J) 

100  CONTINUE 
RETURN 
END 

CCC 

CCC 

SUBROUTINE  SUMF  ( A, B, NBRTIM) 

DIMENSION  A(60,19),  B(19) 

DO  100  J  =  1,19 
B(J)  =  0.0 
DO  100  I  =  1, NBRTIM 
B(J)  =  B(J)  +  A(I,J) 

100  CONTINUE 
RETURN 
END 

CCC 

CCC 

SUBROUTINE  SUMT  (A, IC , NBRTIM) 

DIMENSION  A(60,19),  C(60),  IC(60) 

DO  50  I  =  1,60 
C(I)  =  0.0 
IC(I)  =  0 
50  CONTINUE 

DO  110  I  =  1, NBRTIM 
DO  100  J  =  1,19 

C(I)  =  C(I)  +  A(I,J) 

100  CONTINUE 

IC(I)  =  INT(C(I)) 

110  CONTINUE 
RETURN 
END 
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PROGRAM  CUTDOWN  (INPUT, OUTPUT, TAPEl , TAPES) 

CCC 

CCC  THIS  PROGRAM  READS  IN  TAPEl  (TALKS) 

CCC  WHICH  CONTAINS  64  CHANNELS  OF  DATA  AND 
CCC  COMPRESSES  IT  TO  16  CHANNELS  ACCORDING 
CCC  TO  GUYOTE  AflD  SISSON  1977,  CREATING 
CCC  TAPES  (PERMFILE  -  SIXTEEN) 

CCC 

DIMENSION  A(60,64),  B(60,16) 

NBRTIM  =  60 
DO  100  I  =  1,90 

CALL  READIN  (A,NBRTIM,IWORD , IGRAV, ISPKER, IWRDSZ) 

CALL  CONVERT  ( A, B, NBRTIM) 

CALL  BOUT  (B, NBRTIM, IWORD, IGRAV, ISPKER, IWRDSZ) 

100  CONTINUE 
STOP 
END 

CCC 

CCC 

SUBROUTINE  READIN  (A, NBRTIM, IWORD, IGRAV , ISPKER, IWRDSZ) 

CCC 

CCC  THIS  SUBROUTINE  READS  IN  A  AND  HOLDS 
CCC  THE  BEGINNING  WORD /GRAVITY  MARKER 

CCC 

DIMENSION  A(60,64) 

CCC 

READ  (1,*)  IWORD, IGRAV, ISPKER, IWRDSZ 
DO  100  I  »  1, NBRTIM 
DO  100  J  =  1,64 
READ  (1,*)  A(I,J) 

100  CONTINUE 
RETURN 
END 

CCC 

CCC 

SUBROUTINE  BOUT  (B, NBRTIM, IWORD, IGRAV, ISPKER, IWRDSZ) 

CCC 

CCC  THIS  SUBROUTINE  WRITES  OUT  B  TO  TAPES 
CCC 

DIMENSION  B(60,16) 

CCC 

WRITE  (S,*)  IWORD, IGRAV, ISPKER, IWRDSZ 
DO  100  I  =  1, NBRTIM 
DO  100  J  -  1,16 
WRITE  (S,*)  B(I,J) 

100  CONTINUE 
RETURN 
END 

CCC 

CCC 
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SUBROUTINE  CONVERT(A,B,NBRTIM) 

CCC 

CCC  THIS  SUBROUTINE  CONVERTS  A(60X64)  TO 
CCC  B(60X16)  ACCORDING  TO  THE  LOGRITHMIC 
CCC  COMPRESSION  SCALE  USED  BY  GUYOTE  AND 
CCC  SISSON  1977  AFIT  THESIS 
CCC 

DIMENSION  A(60,64),  B(60,16) 

DO  900  I  =  l.NBRTIM 
DO  100  J  =  1,6 
B(I,J)  =  A(I,J) 

100  CONTINUE 
J  =  7 

DO  210  K  =  7,9 
ADDJ  =  0 
DO  200  N  =  1,2 

ADDJ  =  ADDJ  +  A(I,J) 

J  =  J  +  1 
200  CONTINUE 

B(I,K)  =  ADDJ 
210  CONTINUE 
CCC 

CCC  J  =  13 
CCC 

DO  310  K  =  10,11 
ADDJ  =  0 
DO  300  N  =  1,4 

ADDJ  =  ADDJ  +  A(I,J) 

J  =  J  +  1 
300  CONTINUE 

B(I,K)  =  ADDJ 
310  CONTINUE 
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ccc 

B(I,12)  =  0 
DO  400  J  =  21,25 

B(I,12)  =  B(I,12)  +  A(I,J) 
400  COtJTINUE 
CCC 

B(I,13)  =  0 
DO  500  J  =  26,31 

B(I,13)  =  B(I,13)  +  A(I,J) 
500  CONTINUE 
CCC 

B(I,14)  =  0 
DO  600  J  =  32,40 

B(I,14)  =  B(I,14)  +  A(I,J) 
600  CONTINUE 
CCC 

B(I,15)  =  0 
DO  700  J  =  41,50 

B(I,15)  =  B(I,15)  +  A(I,J) 
700  CONTINUE 
CCC 

B(I,16)  -  0 

DO  800  J  =  51,64 

B(I,16)  =  B(I,16)  +  A(I,J) 
800  CONTINUE 
900  CONTINUE 
RETURN 
END 
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PROGRAM  CONVERT  (INPUT, OUTPUT, TAPE2,TAPE3) 

CCC 

CCC  THIS  PROGRAM  CONVERTS  FROM  THE  64  FREQUENCY 

CCC  FILE  TALKDTA  TO  THE  19  CHANNELS  USED  BY  ^ 

CCC  THRESHOLD  TECHNOLOGY  AS  INPUT  TO  THEIR 
CCC  SPEECH  RECOGNIZER  SYSTEM. 

CCC 

DIMENSION  A(60,64),  B(60,19) 

CCC 

DO  100  K  =  1,90 
NBRTIM  =  60 

CCC 

CALL  READIN  (A,NBRTIM, IWORD , IGRAV, ISPKER, IWRDSZ) 

CALL  NEWFREQ  ( A, B, NBRTIM) 

CALL  WRITEO  (B, NBRTIM, IWORD, IGRAV, ISPKER, IWRDSZ) 

100  CONTINUE 
STOP 
END 

CCC 

CCC 

SUBROUTINE  READIN  (A, NBRTIM, IWORD, IGRAV, ISPKER, IWRDSZ) 
DIMENSION  A(60,64) 

READ  (2,*)  IWORD, IGRAV, ISPKER, IWRDSZ 
DO  100  X  =  1, NBRTIM 

DO  100  J  =  1,64 

READ  (2,*)  A(I,J) 

100  CONTINUE 
RETURN 
END 

CCC 

CCC 

SUBROUTINE  WRITEO  (B, NBRTIM, IWORD, IGRAV, ISPKER, IWRDSZ) 
DIMENSION  B(60,19) 

WRITE  (3,*)  IWORD, IGRAV, ISPKER, IWRDSZ 
DO  100  I  =  1, NBRTIM 

DO  100  J  =  1,19 

WRITE  (3,*)  B(I,J) 

100  CONTINUE 
RETURN 
END 


SUBROUTINE  NEWFREQ  (A,B,NBRTIM) 

CCC 

DIMENSION  A(60,64),  B(60,19) 

CCC 

DO  900  KTIME  =  l.NBRTIM 
DO  900  K  =  1,10 

IF  (K.EQ.l)  FREQ  =  260.0 
IF  (K.EQ.2)  FREQ  =  395.0 
IF  (K.EQ.3)  FREQ  =  535.0 
IF  (K.EQ.4)  FREQ  =  683.0 
IF  (K.EQ.5)  FREQ  =  841.0 
IF  (K.EQ.6)  FREQ  =  1011.0 
IF  (K.EQ.7)  FREQ  =  1198.0 
IF  (K.EQ.8)  FREQ  =  1405.0 
IF  (K.EQ.9)  FREQ  =  1635.0 
IF  (K.EQ.IO)  FREQ  =  1892.0 

IF  (K.EQ.ll)  FREQ  =  2179.0 

IF  (K.EQ.12)  FREQ  =  2505.0 

IF  (K.EQ.13)  FREQ  =  2885.0 

IF  (K.EQ.14)  FREQ  =  3326.0 

IF  (K.EQ.15)  FREQ  =  3855.0 

IF  (K.EQ.16)  FREQ  =  4484.0 

IF  (K.EQ.17)  FREQ  =  5263.0 

IF  (K.EQ.18)  FREQ  =  6277.0 

IF  (K.EQ.19)  FREQ  =  7626.0 

CCC 

B(KTIME,K)  =  0.0 
XO  =«  FREQ  /  4.0 
XI  -  FREQ  /  2.0 
X2  =  FREQ 
X3  =  FREQ  *  2.0 
X4  =•  FREQ  *  4.0 

CCC 

DO  900  J  =  1,64 
X  =  (J  *  78.125) 

FOFX  -  .244  +  ((X-X0)*(1.6888)/FREQ) 

1  -  ((X-X0)*(X-X1)*(1.363)/(FREQ**2)) 

2  +  ((X-X0)*(X-X1)*(X-X2)*(.6265)/(FREQ**3)) 

3  -  ((X-X0)*(X-X1)*(X-X2)*(X-X3)*(.14467)/(FREQ**4)) 

CCC 

IF  (FOFX. GT. 0.0)  B(KTIME,K)  =  B(KTIME,K)  +  A(KTIME, J)*FOFX 

CCC 

900  CONTINUE 
RETURN 
END 
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Explanation  of  Equation  FOFX  =  .  . 


The  frequency  response  for  each  of  Threshold  Technology's  nine¬ 
teen  frequencies  was  the  following  function; 


This  led  to  the  following  table  o'f'  divided  differences: 


X 

f{x) 

^^^n’^n+1) 

^^^n’^n+1 ’^n+2^ 

"Vs’ 

F/4 

. 24444 

1.6888/F 

-1.363/F2 

.6265/F3 

-.144666/F'^ 

F/2 

.6666 

.666668/F 

-.266667/F2  1 

.084/F3 

F 

1.0 

-.266667/F 

-.0277778/F2 

2F 

.7333 

-.183383/F 

1 

4F 

.36666 

1 

1 

Thus,  Newton's  equation  for  f(x)  becomes  the  equation  FOFX  used  in 
subroutine  NEWFREQ. 

(Newton's  formula:  f(x)  =  f[XQ]  +  {x-XQ)f[XQXi  ]  + 

(x-XQ)(x-x^)f[XoX^X2].  .  . 
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PROGRAM  GETDATA  (INPUT, OUTPUT, TAPE9, TAPE 1) 
REWIND  9 
REWIND  1 

CALL  BYPASS  (45) 

CALL  10  (23,30,1,5,3) 

CALL  BYPASS  (41) 

CALL  10  (25,30,1,5,3) 

CALL  BYPASS  (40) 

CALL  10  (26,30,1,5,3) 

CALL  BYPASS  (62) 

CALL  10  (28,30,2,5,3) 

CALL  BYPASS  (42) 

CALL  10  (27,30,2,5,3) 

CALL  BYPASS  (39) 

CALL  10  (30,30,2,5,3) 

CALL  BYPASS  (42) 

CALL  10  (29,30,2,5,3) 

CALL  BYPASS  (56) 

CALL  10  (38,40,3,5,3) 

CALL  BYPASS  (36) 

CALL  10  (36,40,3,5,3) 

CALL  BYPASS  (38) 

CALL  10  (43,40,3,5,3) 

CALL  BYPASS  (39) 

CALL  10  (37,40,3,5,3) 

CALL  BYPASS  (332) 

CCCCCCCCCCC 

CALL  10  (45,40,1,1,2) 

CALL  BYPASS  (81) 

CALL  10  (32,30,2,1,2) 

CALL  BYPASS  (74) 

CALL  10  (46,40,3,1,2) 

CALL  BYPASS  (179) 

CALL  10  (44,40,1,1,2) 

CALL  BYPASS  (78) 

CALL  10  (29,30,2,1,2) 

CALL  BYPASS  (85) 

CALL  10  (44,40,3,1,2) 

CALL  BYPASS  (134) 

CALL  10  (46,40,1,1,2) 

CALL  BYPASS  (80) 

CALL  10  (32,30,2,1,2) 

CALL  BYPASS  (88) 

CALL  10  (39,40,3,1,2) 

CALL  BYPASS  (135) 

CALL  10  (59,40,1,1,2) 

CALL  BYPASS  (64) 

CALL  10  (34,30,2,1,2) 

CALL  BYPASS  (73) 

CALL  10  (40,40,3,1,2) 

CALL  BYPASS  (430) 


ccccccccccc 

CALL  10  (73,60,1,2,2) 
CALL  BYPASS  (37) 

CiVLL  10  (39,40,2,2,2) 
CALL  BYPASS  (50) 

CALL  10  (44,40,3,2,2) 
CALL  BYPASS  (187) 

CALL  10  (60,60,1,2,2) 
CALL  BYPASS  (44) 

CALL  10  (35,40,2,2,2) 
CALL  BYPASS  (48) 

CALL  10  (52,40,3,2,2) 
CALL  BYPASS  (126) 

CALL  10  (48,60,1,2,2) 
CALL  BYPASS  (56) 

CALL  10  (30,40,2,2,2) 
CALL  BYPASS  (64) 

CALL  10  (38,40,3,2,2) 
CALL  BYPASS  (134) 

CALL  10  (52,60,1,2,2) 
CALL  BYPASS  (53) 

CALL  10  (32,40,2,2,2) 
CALL  BYPASS  (62) 

CALL  10  (39,40,3,2,2) 
CALL  BYPASS  (500) 
CCCCCCCCCCC 

CALL  10  (28,30,2,4,2) 
CALL  BYPASS  (49) 

CALL  10  (50,40,3,4,2) 
CALL  BYPASS  (118) 

CALL  10  (51,40,1,4,2) 
CALL  BYPASS  (38) 

CALL  10  (30,30,2,4,2) 
CALL  BYPASS  (47) 

CALL  10  (37,40,3,4,2) 
CALL  BYPASS  (133) 

CALL  10  (41,40,1,4,2) 
CALL  BYPASS  (15) 

CALL  10  (33,30,2,4,2) 
CALL  BYPASS  (37) 

CALL  10  (39,40,3,4,2) 
CALL  BYPASS  (413) 


cccccccccccccccccccccccccccccccccc 

CALL  10  (40,35,4,5,3) 

CALL  BYPASS  (33) 

CALL  10  (35,35,4,5,3) 

CALL  BYPASS  (30) 

CALL  10  (32,35,4,5,3) 

CALL  BYPASS  (34) 

CALL  10  (35,35,4,5,3) 

CALL  BYPASS  (54) 

CALL  10  (20,20,5,5,3) 

CALL  BYPASS  (45) 

CALL  10  (21,20,5,5,3) 

CALL  BYPASS  (45) 

CALL  10  (17,20,5,5,3) 

CALL  BYPASS  (43) 

CALL  10  (18,20,5,5,3) 

CALL  BYPASS  (64) 

CALL  10  (18,20,6,5,3) 

CALL  BYPASS  (49) 

CALL  10  (17,20,6,5,3) 

CALL  BYPASS  (45) 

CALL  10  (18,20,6,5,3) 

CALL  BYPASS  (45) 

CALL  10  (18,20,6,5,3) 

CALL  BYPASS  (69) 

CALL  10  (36,35,7,5,3) 

CALL  BYPASS  (33) 

CALL  10  (37,35,7,5,3) 

CALL  BYPASS  (34) 

CALL  10  (34,35,7,5,3) 

CALL  BYPASS  (31) 

CALL  10  (30,35,7,5,3) 

CALL  BYPASS  (443) 
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ccccccccccccccccc 

CALL  10  (32,30,4,1,1) 
CALL  BYPASS  (112) 

CALL  10  (33,30,4,1,1) 
CALL  BYPASS  (65) 

CALL  10  (30,30,4,1,1) 
CALL  BYPASS  (170) 

CALL  10  (20,20,5,1,1) 
CALL  BYPASS  (76) 

CALL  10  (20,20,5,1,1) 
CALL  BYPASS  (73) 

CALL  10  (20,20,5,1,1) 
CALL  BYPASS  (181) 

CALL  10  (16,20,6,1,1) 
CALL  BYPASS  (86) 

CALL  10  (17,20,6,1,1) 
CALL  BYPASS  (91) 

CALL  10  (16,20,6,1,1) 
CALL  BYPASS  (191) 

CALL  10  (31,30,7,1,1) 
CALL  BYPASS  (65) 

CALL  10  (28,30,7,1,1) 
CALL  BYPASS  (72) 

CALL  10  (28,30,7,1,1) 
CALL  BYPASS  (470) 
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ccccccccccc 

CALL  10  (27,30,4,2,1) 

CALL  BYPASS  (68) 

CALL  10  (28,30,4,2,1) 

CALL  BYPASS  (70) 

CALL  10  (26,30,4,2,1) 

CALL  BYPASS  (73) 

CALL  10  (29,30,4,2,1) 

CALL  BYPASS  (143) 

CALL  10  (21,20,5,2,1) 

CALL  BYPASS  (51) 

CALL  10  (23,20,5,2,1) 

CALL  BYPASS  (61) 

CALL  10  (18,20,5,2,1) 

CALL  BYPASS  (70) 

CALL  10  (18,20,6,2,1) 

CALL  BYPASS  (57) 

CALL  10  (15,20,6,2,1) 

CALL  BYPASS  (58) 

CALL  10  (20,20,6,2,1) 

CALL  BYPASS  (189) 

CALL  10  (24,25,7,2,1) 

CALL  BYPASS  (55) 

CALL  10  (28,25,7,2,1) 

CALL  BYPASS  (51) 

CALL  10  (27,25,7,2,1) 

CALL  BYPASS  (437) 

CALL  10  (17,20,5,3,1) 

CALL  BYPASS  (43) 

CALL  10  (17,20,6,3,1) 

CALL  BYPASS  (35) 

CALL  10  (25,25,7,3,1) 

CALL  BYPASS  (173) 

CALL  10  (20,20,5,3,1) 

CALL  BYPASS  (45) 

CALL  10  (17,20,6,3,1) 

CALL  BYPASS  (99) 

CALL  10  (24,25,7,3,1) 

STOP 

END 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 
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SUBROUTINE  10  (NBREC, IWRDSZ, IWORD, IGRAV, ISPKER) 
DIMENSION  ARRAY  (100,64) 

WRITE  (1,*)  IWORD, IGRAV, ISPKER, ITODSZ 
DO  020  KOUNT  =  1, NBREC 

READ  (9)  (ARRAY(KOUNT, I) ,1=1,64) 

020  CONTINUE 

CALL  MAKEPF  (ARRAY, NBREC, IWRDSZ) 

RETURN 

END 


SUBROUTINE  BYPASS  (NBREC) 
DIMENSION  ARRAY  (64) 

DO  010  KOUNT  =  1, NBREC 

READ  (9)  (ARRAY(I),I=1,64) 
010  CONTINUE 
RETURN 
END 


SUBROUTINE  MAKEPF  (B, NBREC, IWRDSZ) 

THIS  SUBROUTINE  MAKEPF  WILL  CONVERT  ARRAY  B 
(WHICH  CONTAINS  THE  ORIGINAL  UNSCALED  UNMODIFIED 
DATA)  TO  A  STANDARD  SIZE  (IWRDSZ)  AND  THEN  TOITE 
THIS  WORD  BLOCK  TO  A  PERMFILE  FOR  INPUT  TO  SPSS. 

(EACH  WORD  AND  GRAVITY  WILL  HAVE  A  STANDARD  H^^lDSZ; 

ITOTSZ  WILL  BE  THE  MAX  IWRDSZ  AND  WILL  SERVE  AS  A  STANDARD 
TEMPLATE  SIZE  (ZERO  FILLER)  FOR  SPSS  INPUT.  ) 

DIMENSION  8(100,64),  C(100,64) 


ZERO  =0.0 
ITOTSZ  =  60 

IF  (NBREC.lt. IWRDSZ)  CALL 
IF  (NBREC. GT. IWRDSZ)  CALL 
IF  (NBREC.EQ. IWRDSZ)  CALL 
DO  100  I  =  1,  IWRDSZ 

DO  100  J  =  1,64 

WRITE  (1,*)  C(I,J) 

100  CONTINUE 

IF  ( I WRDSZ.GE. ITOTSZ)  GO 
IFILL  =  ITODSZ  +  1 
DO  200  I  =  IFILL,  ITOTSZ 

DO  200  J  =  1,64 

WRITE  (1,*)  ZERO 
200  CONTINUE 
300  RETURN 
END 


ADDTO  (B,NBREC,IWDSZ,C) 
CUT  (B, NBREC, IWRDSZ, C) 
BTOC  ( B, NBREC, C) 


TO  300 
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SUBROUTINE  BTOC  (B,NBREC,C) 

C  THIS  SUBROUTING  COPIES  ARRAY  B  TO  ARRAY  C  WITH  NO 
CHANGES 

DIMENSION  B(100,64),  C(100,64) 

DO  100  I  =  l.NBREC 
DO  90  J  =  1.64 
C(I.J)  =  Bd.J) 

90  CONTINUE 
100  CONTINUE 
RETURN 
END 


SUBROUTINE  ADDTO  (B,NBREC,IWRDSZ,C) 

C  THIS  ROTINE  PADS  ARRAY  C  TO  IWRDSZ  FOR  WRITING 
C  TO  A  PERMFILE  FOR  LATER  INPUT  TO  SPSS. 

C 

DIMENSION  B(100,64),  C(100,64) 

C 

NBRADD  =  IWRDSZ  -  NBREC 

NTH  =  (NBREC  /  (NBRADD  +  1))  +  1 

NHALF  =  NTH  /  2 

IB  =  1 

IC  =  1 

IF  ((NTH*NBRADD).LT. NBREC)  GO  TO  10 
DO  050  J  *  1, NHALF 
DO  45  K  *  1,64 
C(IC,K)  «  B(IB,K) 

45  CONTINUE 
IB  =  IB  +  1 

IC  =  IC  +  I 

50  CONTINUE 
10  CONTINUE 

DO  100  J  =  1,NTH 
DO  90  K  «  1,64 
C(IC,K)  =  B(IB,K) 

90  CONTINUE 
IB  •  IB  +  1 

1C  -  IC  +  1 

IF  (IB. GT. NBREC)  GO  TO  990 
100  CONTINUE 
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C  ADD  AN  NTH  LINE  TO  C  VIA  AN  AVERAGE  OF  B(IB)  AND  B(IB+1) 
C 

C  THERE  WILL  BE  A  LOGICAL  EROR  PROBLEM  BELOW  IF  NBRADD 
C  IS  .GE.  NBREC.  THE  PROBLEM  ISIGNORED  FOR  NOW 
C 

DO  105  K  =  1,64 

r(IC,K)  =  ((B(IB,K)  +  B(IB-1,K))  /  2) 

105  CONTINUE 

IC  =  IC  +  1 
GO  TO  10 
990  CONTINUE 

IF  (IC.GT.IWRDSZ)  GO  TO  999 
DO  110  K  =  1,64 

C  (IC,K)  =  B(NBREC,K) 

110  CONTINUE 
999  CONTINUE 
RETURN 
END 


SUBROUTINE  CUT  ( B, NBREC, IWRDSZ.C) 

THIS  SUBROUTINE  CUTS  DOWN  ARRAY  C  TO  IWRSZ 
FOR  WRITING  TO  PERMFILE  FOR  INPUT  TO  SPSS. 

DIMENSION  B(100,64),  C(100,64) 

NBRCUT  =  NBREC  -  IWRDSZ 
NTH  =  (NBREC  /  (NBRCUT  +1)) 

IB  -  1 
IC  =  1 
10  CONTINUE 

DO  100  J  =  1,NTH 
DO  90  K  =  1,64 
C(IC,K)  =  B(IB,K) 

90  CONTINUE 

IB  =  IB  +  1 
IC  -  IC  +  1 

IF  (IC.GT.IWRDSZ)  GO  TO  999 
100  CONTINUE 

CUT  THE  NTH  LINE  OF  IB 

IB  -  IB  +1 
GO  TO  10 
999  CONTINUE 
RETURN 
END 
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PR0GRA11  TEST( INPUT, OUTPUT, TAPE9,TAPE5=INPUT, 

+  TAPE6=OUTPUT) 

DIMENSION  A(150,68),  FMAT(80),  CMNT(80) 

10  FORMAT(80A1) 

CCC15  F0RMAT(1H1,80A1) 

20  F0RMAT(1H  ,5(A6,F7.2)) 

PRINT  30 

30  F0RMAT(1H  , "INPUT  RTSMAG , IWDGAP, INOPT, IPROPT" , / ) 

READ ( 5 , * )  RTSMAG , IWDGAP , INOPT , IPROPT 
PRINT  35 

35  F0RMAT(1H  , "INPUT  A  FORMAT  OF  UP  TO  80  CHARACTERS",/) 
READ(5,10)  FMAT 
PRINT  40 

40  F0RMAT(1H  , "INPUT  A  TITLE  OF  UP  TO  80  CHARACTERS",/) 
READ(5,10)  CMNT 
IRFRST=1 
IRLAST=0 
ICFRST=1 
ICLAST=64 
IRA=150 
ICA=68 
ICTS=66 
ICSSQ=67 
ICN0RM=68 
INTSRD=0 
INFILE  =  9 
lOFILE  -  6 
WRITE(6,101) 

12  CONTINUE 

CALL  GETTSS( A, IRA, ICA, IRFRST, IRLAST,ICFRST, ICLAST, ICTS , ICSSQ, 
$  INF ILE, FMAT, RTSMAG, IWDGAP, INTSRD) 

IF  (IRLAST.LE. IRFRST)  GOTO  100 
CCCCC  WRITE(6,15)  (CMNT(I) , 1=1 , 80) 

CCCCC  WRITE(6 ,*)  "IRFRST" , IRFRST, "IRLAST" , IRLAST, "ICFRST" , ICFRST, 
CCCCC$  "ICLAST", ICLAST 

CCCCC  WRITE(6,*)  "RTSMAG" , RTSMAG , "IWDGAP" .IWDGAP, "INTSRD" , INTSRD , 
CCCCC$  "INOPT  ".INOPT, "IPROPT", IPROPT 

CALL  NORM{A, IRA, ICA, IRFRST, IRLAST, ICFRST, ICLAST, ICTS , 

$  ICSSQ, ICNORM, INOPT) 

WRITE(6,*)  "SPECTROGRAM" 

CALL  SPECT( A , IRA , ICA, IRFRST , IRLAST .ICFRST , ICLAST , ICTS .ICSSQ , 

$  ICNORM, lOFILE, INOPT, IPROPT) 

GOTO  12 

100  CONTINUE 
WRITE(6,102) 

101  F0RMAT(1HT,"T  SETS  PRINTER  TO  8  LINES  PER  INCH") 

102  F0RMAT(1HS,"S  SETS  PRINTER  TO  6  LINES  PER  INCH") 

STOP 

END 
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FUNCTION  RMAXEL(A, IRA, ICA, IRFRST, IRLAST, ICFRST, ICLAST) 

C****  FUNCTION  RMAXEL(A, IRA, ICA, IRFRST, IRLAST, ICFRST, ICLAST) 

C****  this  function  returns  the  MAXIMUM  ELEMENT  CONTA2NED  IN 
C****  THE  DIMENSION  ARRAY  BETWEEN  ROWS  IRLAST  AND  ICLAST  INCLUSIVE 
C****  BETWEEN  COLUMNS  ICFRST  AND  ICLAST  INCLUSIVE  IN  EACH  ROW 

C****  INPUT 

C****  ARRAY  IS  ANY  DIMENSION 

C****  IRA  IS  THE  TRUE  ROW  SIZE  OF  A 

C****  ICA  IS  AN  APPROPRIATE  COLUMN  SIZE  FOR  A, USUALLY  EQUAL  TO  ICLAST 
C****  IRFRST  IS  THE  STARTIN  ROW  FOR  THE  MAX  ELEMENT  SEARCH 

C****  IRLAST  IS  THE  ENDING  ROW  FOR  THE  MAX  ELEMENT  SEARCH 

C****  ICFRST  IS  THE  STARTING  COLUMN  IN  EACH  ROW  FOR  THE  SEARCH 
C****  ICLAST  IS  THE  ENDING  OLUMN  IN  EACH  ROW  FOR  THE  SEARCH 

DIMENSION  A(IRA,ICA) 

RMAXEL  -  A( IRFRST,!) 

DO  200  I- IRFRST, IRLAST 
DO  100  J-ICFRST, ICLAST 

IF  (A(I,J).  GT  .RMAXEL)  RMAXEL  -  A(I,J) 

100  CONTINUE 

200  CONTINUE 
RETURN 
END 

FUNCTION  SUMEL(ARRAY , IRA, ICA, IRFRST, IRLAST, ICFRST, ICLAST) 

c**** 

C****  FUNCTION  SUMEL(ARRAY, IRA, ICA, IRFRST, IRLAST, ICFRST, ICLAST) 

C****  this  FUNCTION  RETURNS  THE  SUMS  OF  COLUMNS  ICFRST  THROUGH  ICLAST 

C****  OVER  ROWS  IRFRST  THROUGH  IRLAST  INCLUSIVE  IN  DIMENSION  ARRAY 

Q**** 

C****  INPUT 

C****  ARRAY  IS  ANY  DIMENSION 

C****  lAR  IS  THE  TRUE  ROW  SIZE  OF  ARRAY 

C****  lAC  IS  USED  TO  SPECIFY  THE  COLUMN  DIMENSION  OF  ARRAY  AND  SHOULD 
C****  SET  TO  AN  APPROPRIATE  VALUE 

C****  IRFRST  IS  THE  STARTING  ROW  FOR  THE  SUMMATION 

C****  IRLAST  IS  THE  ENDING  ROW  FOR  THE  SUMMATION 

C****  ICFRST  IS  THE  STARTING  COLUMN  IN  EACH  ROW  FOR  THE  SUMMATION 

C****  ICLAST  IS  THE  ENDING  COLUMN  IN  EACH  ROW  FOR  THE  SUMMATION 
DIMENSION  ARRAY<IRA,ICA) 

SUMEL  -  0.0 

DO  200  I-IRFRST, IRLAST 
DO  100  J-ICFRST, ICLAST 

SUMEL  -  SUMEL  +  ARRAY(I,J) 

100  CONTINUE 

200  CONTINUE 
RETURN 
END 
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FUNCTION  SUMSQR(ARRAY , IRA, ICA, IRFRST.IRLAST, ICFRST, ICLAST) 

C****  FUNCTION  SUMSQR(ARRAY,IRA,ICA,IRFRST,IRLAST,ICFRST,ICLAST) 

Qifkifk 

c****  THIS  FUNCTION  COMPUTES  THE  SUM  OF  SQUARES  OF  THE  ELEMENTS 

C****  BETWEEN  COLUMNS  ICFRST  AND  ICLAST  OF  ROWS  IRFRST  THROUGH 

C****  IRLAST  INCLUSIVE  IN  THE  DIMENSION  ARRAY 

C**** 

c****  INPUT 

C****  ARRAY  IS  ANY  DIMENSION 

C****  IRA  IS  THE  TRUE  ROW  DIMENSION  OF  ARRAY 

C****  ICA  IS  THE  APPROPRIATE  COLUMN  DIMENSION  FOR  ARRAY 

C****  IRFRST  IS  THE  STARTING  ROW  FOR  THE  CALCULATION 

C****  IRLAST  IS  THE  ENDING  ROW  FOR  THE  CALCULATION 

C****  ICFRST  IS  THE  STARTING  COLUMN  IN  EACH  ROW  FOR  THE  CALCULATION 

C****  ICLAST  IS  THE  ENDING  COLUMN  IN  EACH  ROW  FOR  THE  CALCULATION 

C**** 

DIMENSION  ARRAY (IRA, ICA) 

SUMSQR  »  0.0 
DO  200  I- IRFRST, IRLAST 
DO  100  J-ICFRST, ICLAST 

SUMSQR  -  SUMSQR  +  ARRAY(I,J)*ARRAY(I,J) 

100  CONTINUE 

200  CONTINUE 
RETURN 
END 
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SUBROUTINE  GETTSS(ARRAY , IRA, ICA, IRFRST, IRLAST, ICFRST, ICLAST, 

$  ICTS , ICSSQ .INFILE , FMAT .RTSMAG , IWDGAP , INTSRD) 

C****  SUBROUTINE  GETTSS(ARRAY, IRA, ICA, IRFRST, IRLAST, ICFRST, ICLAST 
C****  ICTS, ICSSQ, INFILE, FMAT, RTSMAG, IWDGAP, INTSRD) 

(]**** 

c****  this  SUBROUTINE  STARTS  AT  THE  CURRENT  RECORD  OF  THE  INPUT  FILE 
C****  INFILE  AND  READS  A  WORD  INTO  THE  DIMENSION  ARRAY  FROM  INFILE. 
C****  WORD  START  LOCATION  IS  DETERMINED  BY  READING  UNTIL  SOME 
C****  INPUT  TIME  SLICE  HAS  ELEMENTS  WITH  SUM  OF  SQUARES  >  RTSMAG. 
C****  from  this  POINT  THE  WORD  IS  CONSIDERED  TO  CONTINUE  UNTIL 
C****  IWDGAP  CONSECUTIVE  TIME  SLICES  EACH  HAVE  ELEMENTS  WITH  SUM 

C****  OF  SQUARES  <=  RTSMAG.  THE  WORD  END  IS  THEN  SET  AT  TWENTY 

C****  TIME  SLICES  BEFORE  THE  LAST  TIME  SLICE  READ.  STORAGE  IN 
C****  ARRAY  STARTS  AT  ROW  IRFRST.  THE  TIME  SLICE  NL^BER  OF  EACH 

C****  SLICE  WITHIN  INFILE  (BASED  ON  INTSRD)  AND  THE  SUM  OF  SQUARES 

C****  OF  EACH  TIME  SLICE  (THE  ENERGY  OF  THE  TIME  SLICE)  ARE  ALSO 

C****  RETURNED  AS  DESCRIBED  BELOW. 

Qiih-kii 

C****  INPUT 

C****  ARRAY  IS  ANY  DIMENSION 

C****  IRA  IS  THE  TRUE  ROW  DIMENSION  OF  ARRAY 

C****  ICA  IS  A  COLUMN  DIMENSION  FOR  ARRAY  SET  APPROPRIATELY  BUT  MUST 

C****  be  equal  to  (ICLAST-ICFRST)+2  AS  A  MINIMUM  (SEE  OTHER  INPUT 

C****  PARAMETER  DESCRIPTIONS  BELOW) 

C****  IRFRST  IS  THE  ROW  IN  ARRAY  AT  WHICH  STORAGE  IS  TO  BEGIN 
C****  IRLAST  IS  USED  FOR  OUTPUT 

C****  ICFRST  IS  THE  COLUMN  IN  EACH  ROW  OF  ARRAY  AT  WHICH  STORAGE 

C****  OF  COEFFICIENTS  IS  TO  BEGIN 

C****  ICLAST  IS  THE  COLUMN  IN  EACH  ROW  OF  ARRAY  AT  WHICH  STORAGE 
C****  OF  COEFFICIENTS  IS  TO  END  *****NOTE  THAT  NUMBER  OF 

C****  COEFFICIENTS  »  (ICLAST-ICFRST)+1 

C****  ICTS  IS  THE  COLUMN  NUMBER  IN  EACH  ROW  OF  ARRAY  IN  WHICH 
C****  THE  TIME  SLICE  NUMBER  OF  THE  TIME  SLICE  IS  TO  BE  STORED 

C****  ICSSQ  IS  THE  COLUMN  NUMBER  IN  EACH  ROW  OF  ARRAY  IN  WHICH 
C****  THE  SUM  OF  SQUARES  OF  THE  ELEMENTS  OF  THE  TIME  SLICE  IS 

C****  TO  BE  STORED 

C****  INFILE  IS  THE  TAPE  NUMBER  OF  THE  FILE  FROM  WHICH  READING  IS 
C****  TO  TAKE  PLACE 

C****  FMAT  IS  AN  ARRAY  OF  80  CHARACTERS, ONE  CHARACTER  PER  WORD 
C****  WITH  WHICH  TO  READ  THE  INPUT  DATA 

C<r***  RTSMAG  IS  THE  REAL  MAGNITUDE  CUTOFF  FOR  SUM  OF  SQUARES  OF 
C****  ELEMENTS  ON  WHICH  TIME  SLICE  INCLUSION  IN  A  WORD  IS  BASED 

C****  IWDGAP  IS  THE  ESTIMATE  OF  THE  MINIMUM  GAP  (IN  NUMBER  OF 
C****  time  SLICES)  BETWEEN  WORDS 

C****  INTSRD  IS  THE  NUMBER  OF  TIME  SLICES  READ  FROM  INFILE  BEFORE 
C****  this  CALL  TO  GETTSS 


Q^Ht*** 

c**** 

c**** 

Q**** 

Q**** 

Q**** 

Qiiititit 

Q**** 

^**** 

c**** 

^**** 

Q**** 

Q**** 

(^**** 

QAAiNrik 

(]**** 

Q**** 

^**** 

(]**** 

(]**** 

Q**A* 

^**** 

c*** 

c*** 

c*** 

100 

c** 

c** 

c** 

c** 

c** 

c** 

c** 

c** 


OUTPUT 

ARRAY  CONTAINS  A  WORD  READ  FROM  INFILE.  STORAGE  BEGINS  IN  ROW 
IRFRST  AND  ENDS  IN  ROW  IRLAST.  COEFFICIENTS  IN  EACH  ROW 
ARE  STORED  FROM  COLUMN  ICFRST  THROUGH  COLUMN  ICLAST.  THE 
TIME  SLICE  NUMBER  OF  EACH  TIME  SLICE  IS  IN  COLUMN  ICTS  OF 
EACH  ROW  AND  THE  SUM  OF  SQUARES  OF  THE  COEFFICIENTS  IS  IN 
COLUMN  ICSSQ  OF  EACH  ROW. 

IRLAST  IS  THE  ROW  IN  WHICH  STORAGE  OF  THE  WORD  ENDS 

ALL  OTHER  PARAMETERS  ARE  UNCHANGED 

SPECIAL  NOTES 

THE  EXTERNAL  FUNCTION  SUMSQR  IS  USED  FOR  FINDING  THE 

SUM  OF  SQUARES  OF  COEFFICIENTS  IN  EACH  TIME  SLICE 

THE  SIZE  OF  ARRAY  IN  ROWS  SHOULD  BE  AT  LEAST  THE  MAXIMUM 
NUMBER  OF  TIME  SLICES  EXPECTED  IN  A  WORD  PLUS  IWDGAP 
PLUS  IRFRST  MINUS  1 

EACH  TIME  SLICE  IS  ASSUMED  TO  CONTAIN  ICLAST-ICFRST+1 
COEFFICIENTS 

IF  THE  NUMBER  OF  TIME  SLICES  IN  A  WORD  WOULD  CAUSE 

A  SUBSCRIPT  TO  ARRAY  LARGER  THAN  IRA  (ARRAY'S  TRUE  ROW  SIZE) 
THE  ROUTINE  WILL  EXIT  AND  ARRAY  WILL  CONTAIN  ALL  TIME 
SLICES  READ  SO  FAR.  THE  CONDITION  IS  INDICATED  BY 

HAVING  IRLAST=IRA  UPON  SUBROUTINE  EXIT 

IF  AN  EOF  ON  INFILE  IS  REACHED  BEFORE  A  WORD  IS  COMPLETED 
(IN  THE  SENSE  OF  IWDGAP  LIMITS  DESCRIBED)  THE  WORD  IS 
CONSIDERED  TO  HAVE  ENDED  AT  THE  END  OF  THE  FILE  AND 
IRLAST  IS  SET  APPROPRIATELY 

DIMENSION  ARRAY ( IRA, ICA),  FMAT(80) 

THIS  SECTION  READS  FROM  CURRENT  RECORD  TO  BEGINNING  OF  A  WORD 

AS  DETERMINED  BY  SETTING  OF  RTSMAG 

I-IRFRST 

READ(INFILE)  (ARRAY(I,J) ,J=ICFRST, ICLAST) 

CALL  SCALE(ARRAY,IRA,ICA,I) 

THIS  CALL  TO  SCALE  SHOULD  ONLY  BE  MADE 
TO  EQUALIZE  64  CHANNELS  OF  INPUT 
(EQUALIZE  HIGH  AND  LOW  FREQUENCIES) 

BEFORE  THE  SPECTROGRAPH  IS  NORMALIZED 
AND  PRINTED  (LEVINE  -  23  JUN  80) 

IF  (EOF(INFILE) .EQ.l)  GOTO  1000 

INTSRD  -  INTSRD  +  1 

ARRAY(I, ICSSQ)  -  SUMSQR( ARRAY , IRA, ICA, I , I , ICFRST , ICLAST) 

IF  (ARRAY(I, ICSSQ). LT. RTSMAG)  GOTO  100 
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c*** 

C***  FALL  THROUGH  TO  THIS  SECTION  TO  READ  A  WORD 
C***  WHERE  READING  PROCEEDS  TO  WORD  END  PLUS  IWDGAP 
C***  TIMESLICES  WHOSE  MAGNITUDE  IS  LESS  THAN  RTSMAG 

C*** 

C***  FIRST  SAVE  TIME  SLICE  NUMBER  OF  LAST  TIME  SLICE  READ  ABOVE 
C***  SINCE  IT  IS  THE  FIRST  TIME  SLICE  IN  THE  WORD 

C*** 

ARRAY ( I , ICTS ) -FLOAT ( I NTS RD ) 

490  NTSSKP-0 
500  I-I+l 

IF  (I.GT.IRA)  GOTO  1100 

READ(INFILE)  (ARRAY( I , J) , J-ICFRST, ICLAST) 

C** 

CALL  SCALE (ARRAY. IRA, IC A, I) 

C** 

IF  (EOF(INFILE).EQ.l)  GOTO  1000 

INTSRD-INTSRD  +  1 

ARRAY (I, ICTS)  -  FLOAT (INTSRD) 

ARRAY(I,ICSSQ)  -  SUMSQR( ARRAY , IRA, tCA, I , I , ICFRST , ICLAST) 
NTSSKP  -  NTSSKP  +  1 

IF  (ARRAY(I,ICSSQ).GT. RTSMAG)  GOTO  490 
IF  (NTSSKP.lt. IWDGAP)  GOTO  500 
IRLAST-I -IWDGAP 
RETURN 

C*** 

C***  GET  TO  1000  IF  EOF  REACHED  ON  INFILE 
1000  CONTINUE 
IRLAST-I -1 
RETURN 

C*** 

C***  GET  TO  1100  IF  ROW  SIZE  OF  A  WAS  ABOUT  TO  BE  EXCEEDED 
1100  CONTINUE 
IRLAST-IRA 
RETURN 
END 


SUBROUTINE  NORM( ARRAY , IRA, ICA, IRFRST , IRLAST , ICFRST , ICLAST , 

$  ICTS.ICSSQ.ICNORM.IOPT) 

(^**** 

C****  SUBROUTINE  NORM(ARRAY , IRA, ICA, IRFRST, IRLAST, ICFRST, ICLAST, 

C****  ICTS,ICSSQ,ICNORM,IOPT) 

C****  T.ilS  SUBROUTINE  NORMALIZES  A  SECTION  OF  THE  DIMENSION  ARRAY. 

C****  the  SECTION  NORMALIZED  IS  FROM  ROW  IRFRST  THROUGH  ROW  IRLAST 
C****  OVER  COLUMNS  ICFRST  THROUGH  ICLAST.  THE  ROUTINE  EXPECTS  TO 
C****  FIND  THE  SUM  OF  SQUARES  OF  COLUMNS  ICFRST  THROUGH  ICLAST  IN 
C****  COLUMN  ICSSQ  OF  EACH  ROW.  COLUMN  ICTS  OF  EACH  ROW  IS  ASSUMED 
C****  TO  BE  AN  INDENTIFYING  NUMBER  BUT  IS  NOT  USED  IN  THE  CURRENT 
C****  VERSION  OF  THE  SUBROUTINE.  THE  NORMALIZATION  FACTOR  USED  TO 
C****  NORMALIZE  A  PARTICALUR  ROW  OF  A  IS  RETURNED  IN  COLUMN  ICNORM 
C****  OF  EACH  ROW.  THE  NORMALIZATION  PERFORMED  DEPENDS  ON  THE  SETTING  OF 
C****  OF  lOPT  AS  DISCUSSED  UNDER  INPUT  DESCRIPTION 
0**** 

C****  INPUT 

C****  ARRAY  IS  ANY  DIMENSION 

C****  IRA  IS  THE  TRUE  ROW  DIMENSION  OF  ARRAY 

C****  ICA  IS  THE  MAXIMUM  COLUMN  NUMBER  USED  IN  ARRAY  AND  MUST  BE 

C****  EQUAL  TO  AT  LEAST  (ICLAST-ICFRST)+4 

C****  IRFRST  IS  THE  ROW  IN  ARRAY  AT  WHICH  NORMALIZATION  IS  TO  START 

C****  IRLAST  IS  THE  ROW  IN  ARRAY  AT  WHICH  NORMALIZATION  IS  TO  END 

C****  ICFRST  IS  THE  STARTING  COLUMN  IN  EACH  ROW  AT  WHICH  NORMALIZA- 

C****  TION  IS  TO  START 

C****  ICLAST  IS  THE  ENDING  COLUMN  IN  EACH  ROW  AT  WHICH  NORMALIZA- 

c****  tion  is  to  end 

C****  ICTS  IS  A  column  NUMBER  WITHIN  EACH  ROW  ASSUMED  TO  CONTAIN 
C****  an  identifying  NUMBER,  NOT  USED  IN  CURRENT  VERSION  OF 

C****  SUBROUTINE 

C****  ICSSQ  IS  THE  COLUMN  NUMBER  WITHIN  EACH  ROW  WHICH  SHOULD 
C****  CONTAIN  THE  SUM  OF  SQUARES  OF  THE  COLUMNS  ICFRST  THRU  ICLAST 

C****  ICNORM  IS  THE  COLUMN  NUMBER  WITHIN  EACH  ROW  IN  WHli.  '  THE 
C****  NORMALIZATION  FACTOR  USED  FOR  THAT  ROW  IS  TO  BE  STORED 

C****  lOPT  IS  THE  NORMALIZATION  OPTION  AS  FOLLOWS 
C****  2  NORMALIZES  EACH  INDIVIDUAL  ROW  WITHIN  ITSELF 

C****  3  NORMALIZES  OVER  THE  ENTIRE  RANGE  OF  ROWS  IRFRST 

C****  THROUGH  IRLAST 

QAAAA 

C****  OUTPUT 

C****  ARRAY  -  EACH  ROW  OF  ARRAY  CONTAINS  COLUMNS  AS  FOLLOWS: 

C****  NCFRST  THRU  NCLAST  CONTAIN  NORMALIZED  TERMS  (ORIGINAL 
C****  COLUMNS  NORMALIZED  ACCORDING  TO  lOPT  ) 

C****  ICTS  IS  UNCHANGED,  ASSUMED  TO  BE  A  TIME  SLICE  NUMBER 

C****  ICSSQ  CONTAINS  THE  TOTAL  ENERGY  IN  THE  ROW  =  SUM  OF 

C****  RETURNED  ELEMENTS  NCFRST  THRU  NCLAST 

C****  ICNORM  CONTAINS  THE  NORMALIZATION  FACTOR  USED  TO 
C****  NORMALIZE  THAT  ROW 
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c**** 

c****  SPECIAL  NOTES 

C****  VARIABLE  COLUMN  DIMENSION  ICA  MUST  BE  EQUAL  TO  A  MINIMUM  OF 
C****  ((ICLAST-ICFRST)+4)  OR  THE  MAXIMUM  OF(ICLAST,ICFRST,ICTS, 

C****  ICSSQ.ICNORM) ,  WHICHEVER  CAUSES  THE  LARGEST  DIMENSION 

C****  the  column  NUMBERS  OF  ICTS,ICSSQ,  AI4D  ICNORM  SHOULD  ALL  BE 
C****  <  ICFRST  AND/OR  >  ICLAST  (I.E.  SHOULD  NOT  LIE  BETWEEN 

C****  ICFRST  AND  ICLAST  INCLUSIVE) 

C****  IF  THIS  SUBROUTINE  IS  EXECUTED  WITH  lOPT  SET  TO  A  VALUE  OTHER 
C****  THAN  2  OR  3  IT  EXITS  WITHOUT  AFFECTING  ARRAY 

C**** 

DIMENSION  ARRAY (IRA, ICA) 

IF  (IOPT.GT.3)  RETURN 
GOTO  (9999,2000,3000) ,IOPT 

C***  ABOVE  CODE  BRANCHES  TO  APPROPRIATE  OPTION, EXITS  ON  BAD  OPTION 

C*** 

2000  CONTINUE 

G***  IOPT-2, NORMALIZE  EACH  ROW  WITHIN  ITSELF 

DO  2100  I-IRFRST,IRLAST 

ARRAY  ( I ,  ICNORM)=>SQRT(ARRAY  (1 ,  ICSSQ)  ) 

DO  2050  J=ICFRST, ICLAST 

ARRAY(I,J)  =■  ARRAY(I,J)  /  ARRAY( I , ICNORM) 

2050  CONTINUE 

2  ICO  CONTINUE 
RETURN 

C*** 

3000  CONTINUE 

C***  lOPT  -  3  NORMALIZE  OVER  ENTIRE  WORD 

SUMOSUM  =0.0 
DO  3100  I-IRFRST,IRLAST 

SUMOSUM  =  SUMOSUM  +  ARRAY (I, ICSSQ) 

3100  CONTINUE 

R.;ORM  -  SQRT( SUMOSUM) 

DO  3200  I=IRFRST,IRLAST 
DO  3150  J=ICFRST, ICLAST 

ARRAY(I,J)  =  ARRAY(I,J)  /  RNORM 
3150  CONTINUE 

ARRAY (I, ICNORM)  =  RNORM 
3200  CONTINUE 
RETURN 

C*** 

C***  NEXT  IS  COMMON  RETURN/EXIT  FOR  BAD  OPTIONS  WHEN  lOPT 
C****  an  INVALID  OPTION  NUMBER 

9999  RETURN 
END 


SUBROUTINE  SPECT(A, IRA, ICA, IRFRST, IRLAST, ICFRST, ICLAST, 

$  ICTS, ICSSq, ICNORM.IOFILE, INOPT, IPROPT) 

C****  SUBROUTINE  SPECT( A, IRA, ICA, IRFRST, IRLAST, ICFRST, ICLAST, 

C****  ICTS , ICSSQ , ICNORM , lOF I LE , INOPT , IPROPT ) 

DIMENSION  A(IRA,ICA),  SPECLN(66) 

20  FORMATdH  ,F5.0,1X,F10.3,1X,66A1) 

30  F0RMAT(1H+,17X,66A1) 

40  FORMATaH+,35X,l6F6.3) 

C***  ***INITIALIZE  ARRAY  SPECLN, WHICH  HOLDS  THE  SPECTROGRAM  LINES 
NCOEFF  »  ICLAST  -  ICFRST  +  1 
ILNGTH  »  ICLAST  -  ICFRST  +  3 
DO  50  I-l, ILNGTH 
SPECLN(I)  =  "  " 

50  CONTINUE 
SPECLN(l)  = 

SPECLN( ILNGTH)  =  "$" 

C***  ***SET  RNG  OVER  WHICH  COEFFICIENTS  RUN  DEPENDING  ON  THE  NORMALI- 
C***  ***ZATION  OPTION  SPECIFIED,  I.E.  DEPENDING 
C***  ***0N  SETTIN  OF  INOPT 

IF  (INOPT. EQ.l)  RNG  -  100.0 

IF  (IN0PT.EQ.3)  RNG=RMAXEL(A, IRA, ICA, IRFRST, IRLAST, ICFRST, ICLAST) 
C***  ***NOTE  NO  SETTING  FOR  IN0PT=»3,  SETTING  IS  SET  FOR  EACH 
C***  ***INDIVIDUAL  ROW  AT  PRINT  TIME 


c*** 

c***  ***D0  LOOP  THROUGH  1500  LOOPS  THROUGH  ROWS 

C***  ***T0  ACCOMPLISH  CREATION  OF  SPECTROGRAM 

DO  1500  J=IRFRST,IRLAST 
C***  ***PROCESS  FIRST  LINE  OF  SPECTROGRAM 

IF  (INOPT.EQ.2)  RNG=  RMAXEL(A,IRA,ICA, J,J,ICFRST,ICLAST) 

C***  ***SPECTROGRAM  CREATION  NEXT  .FIRST  LINE  OF  TWO 

L=1 

DO  1200  K=ICFRST,ICLAST 
L-L+1 

SPECLN(L)="I" 

IF  (A(J,K).LE.(RNG*0.5))  SPECLN(L)="=" 

IF  (A(J,K).LE.(RNG*0.4))  SPECLN(L)="*" 

IF  (A(J,K).LE.(RNG*0.3))  SPECLN(L)="+" 

IF  (A(J,K).LE.(RNG*0.2))  SPECLN(L)="-” 

IF  (A(J,K).LE.(RNG*0.1))  SPECLN(L)="  " 

1200  CONTINUE 

WRITE(I0FILE,20)  A( J, ICTS) ,A( J, ICSSQ) , ( SPECLN(K) ,K-1 , ILNGTH) 
C***  ***PR0CESS  SECOND  LINE  OF  SPECTROGRAM 

L=1 

DO  1400  K=ICFRST,ICLAST 
L-L+1 

SPECLN(L)="H" 

IF  (A(J,K).LE.(RNG*0.9))  SPECLN(L)="#” 

IF  (A(J,K).LE.(RNG*0.8))  SPECLN(L)="X" 

IF  (A(J,K).LE.(RNG*0.7))  SPECLN(L)="=" 

IF  (A(J,K).LE.(RNG*0.6))  SPECLN(L)="-" 

IF  (A(J,K).LE.(RNG*0.4))  SPECLN(L)-"  " 

1400  CONTINUE 

WRITE( lOFILE , 30)  (SPECLN(K) ,K=1 , ILNGTH) 

C***  ***NOW  PRINT  COEFFICIENTS  AND/OR  OTHER  INFORMATION 
IF  ((IPROPT.EQ.2).  AND  .(NCOEFF.LE.16)) 

$  WRITE(IOFILE,40)  (A(J,K) ,K=ICFRST, ICLAST) 

1500  CONTINUE 
RETURN 
END 

C** 


t 
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SUBROUTINE  SCALE(ARAY,IIRA,I1CA,II) 

C** 

C**  "II"  REPRESENTS  THE  ROW  OF  DATA  JUST  READ  IN 
C**  WHICH  IS  TO  BE  EQUALIZED  TO  THE  SAME  SCALE 
C**  USED  TO  REDUCE  6A  CHANNELS  TO  16 
C**  (LEVINE  23  JUN  80  ) 

C** 

DIMENSION  ARAY(IIRA,IICA) 

C** 

DO  10  K»7,12 

ARAY(II,K)  »  (ARAY(II,K)  -  0)  *  2 
10  CONTINUE 

C** 

DO  20  K-13,20 

ARAY(II,K)  «  (ARAY(II,K)  -  0)  *  4 
20  CONTINUE 

C** 

DO  30  K-21,25 

ARAY(H,K)  »  (ARAY(II,K)  -  0)  *  5 
30  CONTINUE 

C** 

DO  40  K-26,31 

ARAY(II,K)  -  (ARAY(II,K)  -  0)  *  6 
40  CONTINUE 

C** 

DO  50  K-32,40 

ARAY(II,K)  »  <ARAY(ir,K)  -  0)  *  9 
50  CONTINUE 

C** 

DO  60  K-41,50 

ARAY(II,K)  -  (ARAY(H,K)  -  0)  *  10 
60  CONTINUE 

C** 

DO  70  K«51,64 

ARAY(11,K)  -  (ARAY(1I,K)  <•  0)  *  14 
70  CONTINUE 
RETURN 
END 
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